Two-fluid scenario for dark energy models in an FRW universe-revisited 



Bijan Saha^, Hassan Amirhashchi^ , Anirudh Pradhan'^ 

^Laboratory of Information Technologies, Joint Institute for Nuclear Research, 141980 Dubna, Russia 

•^E-mail: bijan@jinr.ru 

■^Young Researchers Club, Mahshahr Branch, Islamic Azad University, Mahshahr, Iran 
E-mail: h. amirhashchi@mahshahriau. ac . ir 

■ ^Department of Mathematics, Hindu Post-graduate College, Zamania-232 331, Ghazipur, India 

E-mail: pradhan.anirudh@gmail.com 

o 

(N 



> 
m 



X 



Abstract 



In this paper we study the evolution of the dark energy parameter within the scope of a spatially homoge- 
^ I neous and isotropic Friedmann-Robertson- Walker (FRW) model filled with barotropic fluid and dark energy 

by revisiting the recent results (Amirhashchi et al. in Chin. Phys. Lett. 28:039801, 2011a). To prevail the 
deterministic solution we select the scale factor a{t) = \/t"e* which generates a time-dependent deceleration 
{^j^ parameter (DP), representing a model which generates a transition of the universe from the early decelerating 

phase to the recent accelerating phase. We consider the two cases of an interacting and non-interacting two- 
fluid (barotropic and dark energy) scenario and obtained general results. The cosmic jerk parameter in our 
derived model is also found to be in good agreement with the recent data of astrophysical observations under 
the suitable condition. The physical aspects of the models and the stability of the corresponding solutions 
are also discussed. 
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1 Introduction 



The use of Type la supernovae as standardized light sources — calibrated candles — led to the observational 
■ discovery of dark energy by two groups in 1998 (Riess et al. 1998, Perlmutter et al. 1999). Before the accelerated 
expansion of the universe was revealed by high red-shift supernovae la (SNe la) observations (Riess et al. 1998, 
Perlmutter et al. 1999), it could hardly be presumed that the main ingredients of the universe are dark sectors. 
The concept of dark energy was proposed for understanding this currently accelerating expansion of the universe, 
and then its existence was confirmed by several high precision observational experiments (Bennett et al. 2003; 
Spergel et al. 2003; Tegmark et al. 2004; Abazajian et al. 2004; Hawkins et al. 2003; Verde et al. 2002), 
especially the Wilkinson Microwave Anisotropy Probe (WMAP) satellite experiment. The WMAP shows that 
dark energy occupies about 73%of the energy of the universe, and dark matter about 23%. The usual baryon 
matter, which can be described by our known particle theory, occupies only about 4% of the total energy of 
the universe. Measurements as of 2008, with the greatest weight coming from the combination of supernovae 
with either cosmic microwave background or baryon acoustic oscillation data, show that dark energy makes up 
72 ±3% of the total energy density of the universe, and its equation of state averaged over the last 7 billion years 
is w = 1.00 ±0.1 (Kowalski et al. 2008). This is consistent with the simplest picture, the cosmological constant, 
but also with a great many scenarios of time varying dark energy or extended gravity theories. In order to explain 
why the cosmic acceleration happens, many theories have been proposed. Although theories of trying to modify 
Einstein equations constitute a big part of these attempt, the mainstream explanation for this problem, however, 
is known as theories of dark energies. It is the most accepted idea that a mysterious dominant component, dark 
energy (DE), with negative pressure, leads to this cosmic acceleration, though its nature and cosmological origin 
still remain enigmatic at present. In the concordance model, the energy content of the Universe is dominated 
by a cosmological constant A ~ 1.7 x 10~^^(ey)'^ such that JIa = A/{3Hq) ~ 0.73. Here Hq denotes the Hubble 
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constant which wc parameterize as Hq = 100 hkm Mpc~^ = 2.1332 h x 10~^'^ey. The second component 
of the concordance model is pressure-less matter with fi^ = Pm/Pc = Pm/i-^HQ/SirG) ~ 0.13//i^, where G is 
Newton's constant. 



The simplest candidate of dark energy is the cosmological constant. It is, however, plagued with the so-called 
coincidence problem and the cosmological constant problem (Weinberg 1989; Sahni and Starobinsky 2000; Pee- 
bles and Ratra 2003; Padmanabhan 2003; Copcland ct al. 2006). Another possible form of dark energy is 
provided by scalar fields. Thus some dynamical scalar field, such as quintessence (Wetterich 1988; Ratra and 
Peebles 1988; Caldwell et al. 1998), phantom (Caldwell 2002; Nojiri and Odintsov 2003, 2005; Sahni and Shtanov 
2003; Xiangyun ct al. 2008), quintom (Elizalde et al. 2004; Feng ct al. 2005; Sctarc 2006; Cai ct al. 2007) 
and K-essence (Alimohammadi and Sadjadi 2006; Wei et al. 2007), are proposed as possible candidate of dark 
energy. However, it is worth mentioned here that for these scalar field models the coincident problem still re- 
mains. Although the two dark components are usually studied binder the assumption that there is no interaction 
between them, one cannot exclude such a possibility. In fact, researches show that a presumed interaction may 
help alleviate the coincident problem (Chimento et al. 2003; Chimento and Pavon 2006). A more comprehensive 
review is provided in Copeland et al. (2006). 

If pd and po are density and pressure, respectively, of the DE can be characterized by the equation-of-state 
(EoE) parameter ojd, defined by uju = ^ which is negative for DE. According to the latest cosmological data 
available, the uncertainties are still too large to discriminate among the three cases w < 1, w = 1, and oj > 1: 
w = 1.04+q;5o (Amanullah et al. 2010). Since the quintessence has the property with the EoS w > — 1 and 
the phantom behaves as lo < —1, wc can speculate that the quintom dark energy is a two-component system 
containing quintessence and phantom. The of DE increases with the increase of scale factor, and both the scale 
factor and the phantom energy density can become infinite at a finite time t, a condition known as the "big 
rip" (Caldwell 2002; Caldwell et al. 2003; Frampton and Takahashi 2003; Ncsseris and Pcrivolaropoulos 2004). 
The closer examination shows that the condition w < — 1 is not sufHcient for a singularity occurrence. First 
of all, a transition phantom cosmology is quite possible. Recently, a new scenario to avoid a future singularity 
has been proposed by Frampton et al. (2011, 2012a, 2012b). In this scenario, the EoS parameter oj < —1, so 
the dark energy density increases with time, but oj approaches —1 asymptotically and sufficiently rapidly that 
a singularity is avoided. This proposed non-singular cosmology was termed as a "little rip" (Brevik et al. 2011; 
Nojiri et al. 2012; Granda and Loaiza 2012; Xi et al. 2012; Astashenok et al. 2012 and references therein). The 
evolution of the little rip cosmology is close to that of ACDM up to the present, and is similarly consistent with 
the observational data. 



The cosmological evolution of a two-field dilation model of dark energy was investigated by Liang et al. 
(2009). The viscous dark tachyon cosmology in interacting and non- interacting cases in non-flat FRW Universe 
was studied by Setare et al. (2009). Recently, Amirhashchi et al. (2011a, 2011b) and Pradhan et al. (2011a) 
have studied an interacting and non-interacting two-fiuid scenario for dark energy models in FRW universe. In 
this report we study the evolution of the dark energy parameter within the framework of a FRW cosmological 
model filled with two fluids (barotropic and dark energy) by revisiting the recent work of Amirhashchi et al. 
(2011a) and obtained more general results. The cosmological implications of this two-fluid scenario will be 
discussed in detail in this paper. In doing so we consider both non-interacting and interacting cases. The 
out line of the paper is as follows: In Sect. 2, the metric and the basic equations are described. Sections 3 
and 4 deal with non-interacting and interacting two-fluid models respectively and their physical significances. 
Physical acceptability and the stability of corresponding solutions are analyzed in Sect. 5. Finally, conclusions 
are summarized in the last Sect. 6. 

2 The metric and basic equations 

We consider the spherically symmetric Priedmann- Robertson- Walker (FRW) metric as 



ds-" = -dr + a^{t) 



1 — fcr^ 



+ r^dfl^ 



(1) 



here a{t) is the scale factor and the curvature constant k is —1, 0, -|-1 respectively for open, flat and closed models 
of the universe. 
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The Einstein's field equations (with SttG = 1 and c = 1) read as 



(2) 



where the symbols have their usual meaning and is the two-fluid energy-momentum tensor consisting of dark 
energy and barotropic fluid. 



In a co-moving coordinate system, Einstein's field equations ([2]) for the line element ([T]) lead to 



Ptot 



and 



Ptot = 3 — :t H — ^ 



(3) 



(4) 



where ptot = Pm +Pd and ptot = Pm + Pd- Here Pm and pm are pressure and energy density of barotropic fluid 
and Pd ^ Pd are pressure and energy density of dark fluid respectively. 

The Bianchi identity — leads to T-^ = which yields 



Ptot + 3-(ptot + Ptot) = 0. 



(5) 



The EoS of the barotropic fluid and dark field are given by 

Pm 



and 



L0„ 



UJD = 



Pm 
PD 

pd' 



(6) 



(7) 



respectively. In the following sections we deal with two cases, (i) non-interacting two-fluid model and (ii) 
interacting two-fluid model. 




k=+l k=-l 



■k=0 



Figure 1: The plot of EoS parameter (lod) Vs. t for non-interacting two-fluid scenario. Here po = Ij = 0.5, 
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3 Non-interacting two-fluid model 



First, we consider that two fluids do not interact with each other. Therefor, the general form of conservation 
Eq. ([SJ leads us to write the conservation equation for the dark and barotropic fluid separately as, 

Pm + 3- (p„ = 0, (8) 
a 

and 

Pi5 + 3-(pd+Pd) = 0. (9) 
a 

Here is, of course, a structural difference between Eqs. (8) and (9). Because Eq. (8) is in the form of ojm which 
is constant and hence Eq. (8) is integrable. But Eq. (9) is a function of wd- Accordingly, pu and pn are also 
function oi ujd- Therefore, we can not integrate Eq. (9) as it is a function of ujd which is an unknown time 
dependent parameter. Integration of Eq. ([5]) leads to 

By using Eq. ([TU)) in Eqs. ([3]) and (|31), we first obtain the pD and pu in term of scale factor a{t) 

PD^^[^ + ^)-Poa-'^'^''-^. (11) 

and ^ 

= - ("2- + ^ + 4 V Poc^ma"^^'+"'"^ (12) 

In literature it is common to use a constant deceleration parameter (Akarsu and Kilinc 2010a, 2010b; Amirhashchi 
et al. 2011c; Pradhan et al. 2011a, 2011b; Kumar and Yadav 2011; Yadav 2011; Kumar and Singh 2011), as 
it duly gives a power law for metric function or corresponding quantity. The motivation to choose such time 
dependent DP is behind the fact that the universe is accelerated expansion at present as observed in recent 
observations of Type la supernova (Riess et al. 1998, 2004; Pcrlmutter et al. 1999; Tonry et al. 2003; Cloc- 
chiatti et al. 2006) and CMB anisotropics (Bennett et al. 2003; de Bernardis et al. 1998; Hanany et al. 2000) 
and decelerated expansion in the past. Also, the transition redshift from deceleration expansion to accelerated 
expansion is about 0.5. Now for a Universe which was decelerating in past and accelerating at the present time, 
the DP must show signature flipping (see the Refs. Padmanabhan and Roychowdhury 2003; Amendola 2003; 
Riess et al. 2001). So, in general, the DP is not a constant but time variable. The motivation to choose the 
following scale factor is that it provides a time-dependent DP. 



Under above motivations we take following ansatz for the scale factor, where increase in term of time evolution 

is 

a(t) = \/Fe*, (13) 

where n is a positive constant. If we put n = 0, Eq. (jl3p reduces to a(t) = -v/e* i.e. exponential law of variation. 
This choice of scale factor yields a non-singular cosmology called a "little rip" which is discussed below. It is 
worth mention here that the solutions in both non-interacting and interacting models do not blow up at any 
given epoch for the choice of this ansatz. It should be remembered that a{t) is a unit less function. In Eq. 
()13|) . a is a function of r = where ti is a constant of unit [Time]. As a result, being a function time, a still 
remains unit less. For simplicity here and further we write a as a function of t with t now being unit less. This 
ansatz generalizes the one proposed in (Amirhashchi et al. 2011a; Pradhan et al. 2012). Recently, the ansatz 
([T3|) is also used by Pradhan & Amirhashchi (2011a) in studying the accelerating DE models in Bianchi type-V 
space-time. 



We define the deceleration parameter q as usual, i.e. 



9 = -— = TH- 14 



Using Eq. ^ into Eq. we find 

q=T^^-l. (15) 
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Figure 3: The plot of deceleration parameter (q) Vs. t 
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From Eq. (IT5|) . we observe that q > for t < v2n — n and g < for t > y/2n — n. It is observed that for 
< n < 2, our model is evolving from deceleration phase to acceleration phase. Also, recent observations of 
SNe la, expose that the present universe is accelerating and the value of DP lies to some place in the range 
— 1 < g < 0. It follows that in our derived model, one can choose the value of DP consistent with the observation. 
Figure 3 depicts the deceleration parameter (g) versus time which gives the behaviour of q from decelerating to 
accelerating phase for different values of n. 



By using this scale factor in Eqs. (fTTj) and ([T2|). the pjj and pd are obtained as 



and 



Pd 



n + 1 
2t 



n 



k 



poWm(i"e*)~ 



(l+c^™) 



(16) 



(17) 



respectively. By using Eqs. (|16p and (|17|) in Eq. ([7|), we can find the equation of state of dark field in term of 
time as 

o ( n±t\^ n I k ^ /'J-rl„^^-|(l+LJ,„) 



UJd = 



^ 2t J 



k 



n-f (i+i^,„) 



(18) 



The behavior of EoS for DE {ojd) in term of cosmic time t is shown in Fig. 1. It is observed that for closed 
universe the ojd is decreasing function of time whereas for open and fiat universes the EoS parameter is an 
increasing function of time, the rapidity of their growth at the early stage depends on the type of the universes, 
while later on they all tend to the same constant value independent to it. From this figure we also observe 
that the EoS parameters of closed, open and flat universes are varying in quintessence {ujd > ^l)j phantom 
(—3 < COD < —1) and Super phantom {ujo < —3) regions respectively, while later on they tend to the same 
constant value —1 (i.e. cosmological constant) independent to it. 

From Fig. 1, it is observed that in the case of open and fiat universes w^i is less than —1, so the DE density 
increases with time, but ujo approaches —1 asymptotically and sufficiently rapidly that a singularity is avoided. 
This is the scenario of a "httle rip" Frampton et al. (2011, 2012a, 2012b). 




■ k=+l k=-l 



k=0 



Figure 4: The plot of EoS parameter Vs. t for interacting two-fluid scenario. Here po = 1, Wm = 0.5, n = ^, 
a = 0.3 
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Figure 6: The plot of sound speed Vg Vs. t for interacting two-fluid scenario 
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Figure 7: The plot of energy conditions versus t in non-interacting two-fluid model 

The expressions for the matter-energy density n„i and dark-energy density fl/j are given by 

Prn _ 4pot^ 
3iJ2 ~ 3(t + 



^^n^^ = ,r^^{t-er-^^+--\ (19) 



and 

respectively. Eqs. p9| and ([20)1 reduce to 

From the right hand side of Eq. (PT|) it is clear that in flat universe {k — 0), ft — I and in open universe 
{k = —1), n < 1 and in closed universe (A; = -1-1), £7 > 1. But at late time we see for all flat, open and closed 
universes i7 — > 1. This result is also compatible with the observational results. Since our model predicts a flat 
universe for large times and the present-day universe is very close to flat, so the derived model is also compati- 
ble with the observational results. The variation of density parameter with cosmic time has been shown in Fig. 2. 

In the Universe nearly 70% of the energy is in the form of dark energy. Baryonic matter amounts to only 
3 — 4%, while the rest of the matter (27% is believed to be in the form of a non-luminous component of non- 
baryonic nature with a dust- like equation of state (w = 0) known as cold dark matter (CDM). In this case, 
if the dark energy is composed just by a cosmological constant, then this scenario is called A-CDM model. A 
convenient method to describe models close to A CDM is based on the cosmic jerk parameter j, a dimensionless 
third derivative of the scale factor with respect to the cosmic time (Chiba and Nakamura 1998; Sahni 2002; 
Blandford et al. 2004; Visser 2004, 2005). A deceleration-to-acceleration transition occurs for models with a 
positive value of jo and negative go- Flat A CDM models have a constant jerk j = 1. The jerk parameter in 
cosmology is defined as the dimensionless third derivative of the scale factor with respect to cosmic time 

m = , (22) 



and in terms of the scale factor to cosmic time 
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Figure 8: The plot of energy conditions versus t in non-interacting two-fluid model 



where the 'dots' and 'primes' denote derivatives with respect to cosmic time and scale factor, respectively. The 
jerk parameter appears in the fourth term of a Taylor expansion of the scale factor around ag 



^ = 1 + i/o(t - to) - \qoHl[t - tof + IjoH^it - tof + O [{t - h)^] , 
ao 2 



where the subscript shows the present value. One can rewrite Eq. ((22|) as 



Eqs. (dni) and (gl]) reduce to 



Jit) 



1 



6n 



8n 



(24) 

(25) 
(26) 



This value is overlap with the value j ~ 2.16 obtained from the combination of three kinematical data sets: the 
gold sample of type la supernovae (Riess et al. 2004), the SNIa data from the SNLS project (Astier et al. 2006), 
and the X-ray galaxy cluster distance measurements (Rapetti et al. 2007) for 



50n 



(27) 



where 



A = 0.03 (84100n + 1450 Vl450n3 + 3364^2^ ' . 

4 Interacting two-fluid model 

Secondly, we consider the interaction between dark energy and barotropic fluids. For this purpose we can write 
The continuity equations for dark fluid and barotropic fluids as 



Pni + 3-(Pm, +Pm) = Q, 

a 



and 



Pd + 3-{pD +Pd) 
a 



-Q- 



(28) 
(29) 



The quantity Q expresses the interaction between the dark energy components. Since we are interested in an 
energy transfer from the dark energy to dark matter, we consider Q > 0. Q > 0, ensures that the second law 
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Figure 9: The plot of energy conditions versus t in non-interacting two-fluid model 



of thermodynamics is fulfilled (Pavon and Wang 2009). Here we emphasize that the continuity Eqs. (|28p and 
()29|) imply that the interaction term (Q) should be proportional to a quantity with units of inverse of time i.e 
Q (X J. Therefore, a first and natural candidate can be the Hubble factor H multiplied with the energy density. 
Following Amendola et al. (2007) and Gou et al. (2007), we consider 



Q = 3Hapm, 

where cr is a coupling constant. Using Eq. (15111) in Eq. (^5)) and after integrating, we obtain 



By using Eq. ([31]) in Eqs. ([3]) and (|4]), we again obtain the po and pn in term of scale factor a(t). 



PD 



poa 



-3(H 



and 



PD 



a a 



k 



Po(^™-fT)a-3(i+-'"-), 



a a a 

respectively. Putting the value of a{t) from Eq. (fT3|) in Eqs. (|32|) and ((33)) . we obtain 

2 

I n + T, \ 

Pd = 3 



2t 



3k 



and 



PD 



n + t 
2t 



-^ + ^-Po(^™-a)(tV)-i(^+--'^) 



respectively. Using Eqs. ([M|) and (|55|) in Eq. (O, we can find the EoS parameter of dark field as 



3 (^)' - f + - Poi^m - a)(t»e*)-i(i+"^-'^) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



The behavior of EoS for DE {ojd) in term of cosmic time t is shown in Fig. 4. It is observed that for closed 
universe the ujo is decreasing function of time whereas for open and flat universes the EoS parameter is an 
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Figure 10: The plot of energy conditions versus t in interacting two-fluid model 



increasing function of time, the rapidity of their growth at the early stage depends on the type the universes, 
while later on they all three tend to the same constant value independent to it. From this figure we also observe 
that the EoS parameters of closed, open and flat universes are varying in quintessence {lod > ~1), phantom 
(—3 < LU]j < —1) and Super phantom (lod < —3) regions respectively, while later on they tend to the same con- 
stant value —1 (i.e. cosmological constant) independent to it. From Fig. 4, we also observe that the interaction 
pushes all closed, open and flat universes to darker regions. In this case we also observe that there is a "little 
rip" in the evolution of the open and flat universes. 

The expressions for the matter-energy density and dark-energy density u are given by 



and 



nn = l + ^ - -i^(re*)-i(i+"'"-), (38) 



respectively. From Eqs. (f37| and (pSj) . we obtain 



4fc 

"-"'" + "^ = ^+ t»-V(n-^0- - ^''^ 

which is the same as Eq. (j2ip . Therefore, we observe that in interacting case the density parameter has the 
same properties an non-interacting case. The expressions for deceleration parameter and jerk parameter are also 
same as in the case of non-interacting case. 

Studying the interaction between the dark energy and ordinary matter will open a possibility of detecting 
the dark energy. It should be pointed out that evidence was recently provided by the Abell Cluster A586 in 
support of the interaction between dark energy and dark matter (Bertolami et al. 2007; Delliou et al. 2007). 
We observe that in non-interacting case both open and flat universes can cross the phantom region whereas in 
interacting case only open universe can cross phantom region. 



5 Physical acceptability and stability of solutions 

For the stability of corresponding solutions in both non-interacting and interacting models, we should check that 
our models are physically acceptable. For this, firstly it is required that the velocity of sound should be less 
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Figure 11: The plot of energy conditions versus t in interacting two-fluid model 



than velocity of light i.e. within the range < Vg — ^ 1- 

In our non-interacting and interacting models, we obtained the sound speeds as 



3^ 0.5-l-t ^ I 3 (Q.5+tf 1 I 1 k , _fe l-125(^^ + v^e') 

2i^j + 2^^ F + 2x^ + yte' i:^^ ^^^^ 



and 



3^ 0.5+t \ _ 3 {a.5+ty _ 3 k 3fe_ I 9£Wj2^ 



_3/0^N , 3 (0^5+t£ _ 1 1 fc _fe 0.36(1 ^-H%/te') 

" ~ 1 , /T t ' y*'^) 

3/ 0.5+t \ 3 3 fc __3^ , 9l^j2_ 

respectively. In both cases we observe that Vs < I. From Figures 5 & 6, we observe that in both non- interacting 
and interacting cases Vg < 1. 

Secondly, the weak energy conditions (WEC) and dominant energy conditions (DEC) are given by 

(i) Peff > 0, (ii) Peff - Peff > and (in) p^ff + Peff > 0. 

The strong energy conditions (SEC) are given by p^/f + ipeff > 0. 

From the Figures 7 — 12, we observe that 

• The WEC and DEC for the closed universe in both non-interacting and interacting cases are satisfied. 

• In both open & flat models, the WEC and DEC are satisfied in initial stages of the evolution of the universe 
(i.e. in decelerating phase). In these both models WEC and SEC are violated at later times as expected 
but DEC does not violet. 

• The SEC for both non-interacting and interacting cases are satisfied in early stages of the evolution of the 
universe whereas it violet at present epoch due to acceleration for all three open, closed and flat models 
as expected. 
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Figure 12: The plot of energy conditions versus t in interacting two-fluid model 

Therefore, on the basis of above discussions and analysis, our corresponding solutions are physically acceptable. 

A rigorous analysis on the stability of the corresponding solutions can be done by invoking a perturbative 
approach. Perturbations of the fields of a gravitational system against the background evolutionary solution 
should be checked to ensure the stability of the exact or approximated background solution (Chen and Kao 
2001). Now we will study the stability of the background solution with respect to perturbations of the metric. 
Perturbations will be considered for all three expansion factors Oj via 



UBi + Sui = aBi{l + Sbi) 



(42) 



We will focus on the variables 6bi instead of Sa.i from now on for convenience. Therefore, the perturbations 
of the volume scale factor Vb ~ li^^iUi, directional Hubble factors 9i — ^ and the mean Hubble factor 



Si=3 '3 ~ W '^^^ shown to be 

V ^VB + VBY,Sb,, 



One can show that the metric perturbations Sbi, to the linear order in 

i 

Vb ■ ■s-^ ■ 

■ -^5h + y^5b,eB^^Q, 
Vb ^ 



obey the following equations 



^<55, =0. 



From above three equations, we can easily find 



Sk + 



0, 



(43) 

(44) 

(45) 
(46) 

(47) 



where Vb is the background volume scale factor. In our case, Vb is given by 



Vb 



(48) 
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Figure 13: The plot of the "actual" fluctuations Sui versus t 

Using above equation in equation (47) and after integration we get 

3 3 /3 3 13\ 

= Cit"3"e"4'WittakerM f --n,--n-F -, j , (49) 

where Cj is an integration constant. Therefore, the "actual" fluctuations for each expansion factor 5ai = asiSbi 
is given by 

(3 3 1 3 \ 
"i'^' "4" 2' 2 / ' ^^^^ 

From above equation we see that for n >> 1, Sai approaches zero. Fig. 13 is the plot of the "actual" fluctuations 
Sui versus t which also shows that (5ai as t 00. Consequently, the background solution is stable against 
the perturbation of the graviton field. 

6 Concluding remarks 

In this present work we continue and extend the previous work of Amirhashchi et al. (2011a). In summary, we 

have studied a system of two fluid within the scope of a spatially homogeneous and isotropic FRW model. The 
role of two fluid either minimally or directly coupled in the evolution of the dark energy parameter has been 
investigated. The scale factor is considered to be a power law function of time which yields a time dependent 
deceleration parameter. It is observed that in an interacting and non-interacting cases both open and flat uni- 
verses can cross the phantom region. It is observed that the closed universe is corresponding to quintessence 
whereas the flat and open universes are corresponding to phantom model of universe. During the evolution of 
the universe, we find that the EoS parameter for closed universe changes from ii;>— ltoii;<— 1, which is 
consistent with recent observations. If we put n = 1 in the present paper, we obtain all results of recent paper 
of Amirhashchi et al. (2011a). 

It is observed that the EoS parameters of closed, open and flat universes are varying in quintessence 
{lod > —1), phantom (—3 < ojd < —1) and Super phantom {ojn < —3) regions respectively, while later on they 
tend to the same constant value —1 (i.e. cosmological constant region) independent to it. 

Our special choice of scale factor yields a time dependent deceleration parameter which represents a model of 

Universe which takes evolution from decelerating to accelerating phase which is in good agreement with current 
observations. It is worth mentioned here that for different choice of n, we can generate a class of DE models in 
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FRW universe. It is also observed that such DE models are also in good harmony with current observations. 

We also observe that our corresponding solutions are physically acceptable and the solutions are stable. 
Thus, the sohitions demonstrated in this paper may be useful for better understanding of the characteristic of 
DE in the evolution of universe within the framework of FRW. 
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